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Abstract. We study the thermal quantum correlations and entanglement in spin- 
1 Bose-Hubbard model with two and three particles. While we use negativity to 
calculate entanglement, more general non-classical correlations are quantified using 
a new measure based on a necessary and sufficient condition for zero-discord state. 
£h ' We demonstrate that the energy level crossings in the ground state of the system are 

signalled by both the behavior of thermal quantum correlations and entanglement. 

7— I 1 

" 1. Introduction 

Entanglement, being considered as the resource of quantum information science, has 
been utilized to investigate various properties of condensed matter systems [1,2]. 
However, it has been discovered that entanglement is not the only kind of useful non- 
classical correlation present in quantum systems. Numerous quantifiers of quantum 
correlations have been proposed to reveal the non-classical correlations that cannot be 
captured by entanglement measures [3]. 

Quantum phase transitions (QPTs) are critical changes in the ground states of 
many-body systems when one or more of its physical parameters are continuously 
changed at absolute zero temperature. These qualitative changes can be signalled by 
energy level crossings in the ground state of the system. Signatures of QPTs are also 
present at sufficiently low but experimentally accessible temperatures, where thermal 
fluctuations are still negligible as compared to quantum fluctuations. In recent years, the 
techniques of quantum information theory have been implemented for the investigation 
of quantum critical systems. For instance, several correlation measures have been shown 
to identify the critical points of QPTs with success in quantum critical spin chains [3-5]. 

In this work, we analyze quantum correlations in a spin-1 Bose-Hubbard model 
with two and three particles considering periodic boundary conditions. As a measure of 
quantum correlations, we use a recently introduced measure for an arbitrary bipartite 
system based on a necessary and sufficient condition for a zero-discord state in the 
coherence- vector representation of density matrices [6]. On the other hand, we adopt 
negativity to measure the amount of entanglement in a quantum state. We demonstrate 
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that both quantum correlations that are more general than entanglement and negativity 
mark the critical points corresponding to energy level crossings in the ground state of 
the system. Although we only consider few particle systems in our study, this interesting 
behavior might have consequences even for actual quantum critical systems, where the 
number of particles is very large and the energy level crossings really lead to quantum 
QPTs. 

2. Correlation Measures 

In this section, we briefly review the measures of quantum correlations used in our work. 
We first introduce a measure of non-classical correlations recently proposed by Zhou et 
al. based on a necessary and sufficient condition for a zero-discord state [6]. A general 
bipartite state p ab can be expressed in coherence-vector representation as 

ab = _^ja ^ jb + J2 Xi Xi ® L + <g) J2 VjYj 

mn ^-^ 2n 2m ^-^ 

1=1 3=1 

^ m 2 -l n 2 -l 
i=l j=l 

where the matrices {Xi : % — 0, 1, • • • , m 2 — 1} and {Yj : j = 0, 1, • • • , n 2 — 1}, satisfying 
tr(X k Xi) = tr(Y k Yi) = 25 k i, define an orthonormal Hermitian operator basis associated 
to the subsystems a and b, respectively. Here, / is the identity matrix for the specified 
subsystem. The components of the local Bloch vectors x = {xi}, y = {yj} and the 
correlation matrix T = tij can be obtained as 

Xi =trp a \X t ®I b ), 

Vj =trp a6 (r»®^-), 

U^tTp^Xi^Yj). (2) 

By making use of the above representation of bipartite quantum states, the measure of 
non-classical correlations is given by 

m 2 -l 

Q(p ab ) = 4 E ( 3 ) 

i=m 

where A, are the eigenvalues of the criterion matrix A = TT l — y^yxx 1 in decreasing 
order. The motivation behind the definition of this measure and details of its derivation 
can be found in Ref. [6]. 

Negativity is a measure of entanglement that can be straightforwardly calculated for 
an arbitrary bipartite system in all dimensions. Although we cannot conclude whether 
a positive partial transpose state (zero negativity state) is entangled or separable in 
general, negativity is still a reliable measure for all negative partial transpose states [7,8]. 
For a given bipartite density matrix p ab , it can be defined as the absolute sum of the 
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negative eigenvalues of partial transpose of p ah with respect to the smaller dimensional 

system, 

N(p ab ) = \j2^\-^ (4) 

i 

where rji are all of the eigenvalues of the partially transposed density matrix (p ab ) tA . 



3. Spin-1 Bose-Hubbard Model 

The Hamiltonian describing the system of spin-1 atoms in an optical lattice is given by 
[9,10] 

+ yE(( S ^) 2 -2n l ), (5) 

i 

where a ia (a^) is the creation (annihilation) operator for an atom on site i with z 
component of its spin being equal to o = —1,0,1. Here hi = Y2a a l<T a iv * s the total 
number of atoms on site % and S l tot gives the total spin on ith lattice site. The parameter 
t represents the tunneling amplitude, Uq is the on-site repulsion and U2 differentiates 
the scattering channels between atoms with different S to t values. 

From this point on, we assume that the temperature is low enough and the tunneling 
amplitude t is small so that the overlap between the wavefunctions of the particles in 
neighboring sites is almost zero. Under these assumptions, the spin-1 Bose-Hubbard 
Hamiltonian can be treated perturbatively. Second order perturbation theory in t gives 
the effective Hamiltonian as [10] 

In addition to the original spin-1 Bose-Hubbard Hamiltonian, an external magnetic field 
u has been added to the effective Hamiltonian. Sj is the spin operator of the particle 
on site % with J = ^ Si and / represents the identity operator. In terms of the original 
Bose-Hubbard Hamiltonian parameters t, U , U 2 , the effective coupling constants r, r, 
7 for single particle per site are given by 

At At 

r ~ Wo + U 2 ) ~ 3{U -2U 2 y 
2t 



T = 



C/0 + C/2' 



= 2t + (7) 

with r = r — 7. In what follows, we will consider the two and three particle cases with 
a single particle per site. 
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3.1. Two particles 

In this setting, the explicit form of the effective Hamiltonian given by Eq. (6) reads 

# 2 = wJ z + r/ + TS 1 -S 2 + 7 (Si-S 2 ) 2 . (8) 

Using the identity Si • S 2 = (J 2 — Sf — Sf)/2, the two particle Hamiltonian H 2 can be 
written in the total spin basis as 

H 2 = co.J z + \{J 2 ~ 4/) + \{J 2 - 4/) 2 + rl, (9) 

where the energy eigenvalues are determined as E JM = wM+r(j(j + l)-2)/2 + 7[(j(j + 
1) — 4) 2 — 4]/4. The density matrix of our system at finite temperature T can be written 

as 

e -/3H 

PT = (10) 

with the partition function of the system is given by Z = tr(e~" H ) = e~^ T [2 cosh /3t(1 + 
2 cosh f3u) + e -Pi^- 2T ) + 2e _/3r cosh 2/3r] and (3 = 1/T with Boltzmann constant k B = l. 

In Fig. 1 (a) and (b), we present our results related to the thermal entanglement and 
quantum correlations in the system of two particles as a function of r when 7 = ou — 1 
for T = 0.05, 0.5, 1. Leggio et al. have recently discussed the behavior of thermal 
entanglement in this model, revealing a connection between the different phases of 
entanglement and the energy level crossings in the ground state of the system [11]. 
We demonstrate here that not only the negativity but also the non-classical correlations 
of the system experience two sharp transitions at points r = 0.5 and r = 4 when the 
temperature is sufficiently low. Examining the Fig. 1 (c), it is not difficult to see that 
these sharp transitions are connected with the appearance of energy level crossings in 
the ground state of the system. In fact, ground state crossings occur at the points 
2r = u) and r = u> + 3j, and the connection between the crossings and the considered 
correlation measures is independent of the values of 7 and oj. We also note that when 
r < 0, non-classical correlations in the system grows and reaches to a constant value in 
this regime as the temperature is increased. 

3.2. Three particles 

When it comes three particles, the effective Hamiltonian with periodic boundary 
conditions takes the form 

H 3 = uoJ z + rl + r(Si • S 2 + S 2 • S 3 + S 3 • Si) 

+ 7[(Si-S 2 ) 2 + (S 2 -S 3 ) 2 + (S3-Si) 2 ]. (11) 

Similarly to the case of two particles, we straightforwardly obtain the energy eigenvalues 
of the Hamiltonian and the thermal density matrix p? to evaluate the negativity 
and non-classical correlations in the system. In this case, negativity and quantum 
correlations are calculated considering the bipartition of 3 ® 9, that is, we look at the 
correlations between the first particle and the remaining two particles in the system. 
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Despite the fact that we do not investigate the multipartite non-classical correlations, 
one can indeed use tripartite negativity defined in Ref. [12] to analyze the multipartite 
entanglement. It is easy to see that, due to the symmetry of the considered system, 
the tripartite negativity reduces to usual negativity which is calculated by taking the 
partial transpose with respect to any of the three qubits. 

Fig. 2 (a) and (b) display our results on the thermal entanglement and quantum 
correlations in the system of three particles with periodic boundary conditions as a 
function of r when 7 = uj = 1 for T = 0.05, 0.5, 1. The low lying energy levels and their 
crossing points are also shown in Fig. 3 (c). Looking at the figures, we observe that the 
two sudden jumps of negativity and quantum correlations correspond to the crossings 
of the energy levels in the ground state of the system at r = 1/3 and r = 2/3. We 
note that, different from the case of two particles, negativity and quantum correlations 
do not show a decreasing behavior about the second transition point, r = 2/3, in case 
of three particles. Moreover, the plateau occurring after the first transition here is 
considerably shorter as compared to the two particle case. Lastly, the revival of non- 
classical correlations with increasing temperature can also be seen when r < 0. 

4. Conclusion 

In summary, we have investigated the thermal quantum correlations and entanglement 
in a spin-1 Bose- Hubbard model with two and three particles with periodic boundary 
conditions. Our results demonstrate that both the behavior of thermal quantum 
correlations and entanglement spotlight the energy level crossings in the ground 
state of the system. Despite the fact that our discussion is limited to few particle 
systems, the connection between the behavior of correlations measures and ground 
state crossings might have consequences even for real quantum critical systems having 
large number of particles. Finally, it would be interesting to analyze the relation of 
some thermodynamical quantities (such as the specific heat) and correlations measures 
since various non-trivial behaviors of certain thermodynamical quantities might give 
information about the correlations in the system [13,14]. 
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Figure 1. The thermal entanglement (a) and quantum correlations (b) of Spin-1 Bose- 
Hubbard model with two particles as a function of the parameter r when 7 = uj = 1 
for T = 1 (dotted line), T = 0.5 (dashed line) and T = 0.05 (solid line). The low lying 
energy levels and their crossings in the ground state of the system are displayed in (c) . 



8 



1.0 
0.8 

*-t— » 

US 

jj 0.4 

z 

0.2 
0.0 



(a) 



-4 



-2 



/ ' 







| 0.6 : (b) 

| o.5 ; 
go.4; 

"0.3 1 
s 

B 0.2 : 

o i 



o.o 



-4 



-2 




r 



(50 

a 
w 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 



Figure 2. The thermal entanglement (a) and quantum correlations (b) of Spin-1 Bose- 
Hubbard model with three particles as a function of the parameter r when 7 = uj = 1 
for T = 1 (dotted line), T = 0.5 (dashed line) and T = 0.05 (solid line). The low lying 
energy levels and their crossings in the ground state of the system are displayed in (c) . 



